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ABSTRACT 



Wc investigate the orbital dynamics of hierarchical three-body systems containing 
a double neutron star system orbiting around a massive black hole. These systems 
show complex dynamical behaviour because of relativistic coupling between orbits of 
the neutron stars in the double neutron star system and the orbit of the double neutron 
star system around the black hole. The orbital motion of the neutron stars around 
each other drives a loop mass current, which gives rise to gravito-magnetism. Generally, 
gravito-magnetism involves a rotating black hole. The hierarchical three-body system 
that we consider is an unusual situation in which black hole rotation is not required. 
Using a gravito-electromagnetic formulation, we calculate the orbital precession and 
nutation of the double neutron star system. These precession and nutation effects arc 
observable, thus providing probes to the spacetimc around black holes as well as tests 
of gravito-electromagnetism in the framework of general relativity. 

Key words: black hole physics - gravitation - gravitational waves - pulsars: general 
- stars: neutron binaries: close 



1 INTRODUCTION 



. ' There is strong evide nce that the Galac tic Centre contains a black hole with mass M « 4 x 10 6 M Q (|Gillessen et alj|2009ft . X- 
\ ray observations (e.g. iMuno et al.ll2009T ) indicate that the Galactic central region contain s a large number of stellar remnants. 
?H ■ There have been studies su ggesting that about 20,000—40,000 stella r-mass black holes (|Miralda-Escude fc Gouldll2000l ) and 
thousands of neutron stars (|Freitag et aT1l2006l ; IWharton etai1l2012l ) are residing in the region. The pr esence of such a large 
black hole population requires tha t the Galactic Centre is dynamically relaxed. However, recent studies (|Buchholz et al.l l 2009: 
IDo et al.ll2009l : [Bartko et alj|201(if) have ind i cated otherwise, implyin g that there may well be much fewer stellar-mass black 



holes in the Galactic Centre (|Merrittj[2uTol : lAntonini fc Peretsll20r2T ). The estimate of ~ 10 3 neutron stars is based on the 



models assuming the presence of a cluster of ~ 10 stellar-mas s black holes in th e same region. Without this stellar-mass 



black hole cluster, the number of neutron stars would be higher (Frcitag et al. 2006). A half dozen pulsars have already been 
discovered in the central parsec of Sgr A* (|Deneva et al.l [20091 : iMacquart et alj|2010h . showing the evidence of the neutron 
star population in the Galactic Centre. Centres of galaxies with spheroids similar to that of our Galaxy are expected to host 
thousands or more neutron stars. Large elliptical galaxies have massive stellar spheroids and their neutron star populations 
are scaled accordingly. Some of these neutron stars would eventually fall into the massive blac k holes, forming EMRI (extrem e 
mass-ratio inspiral) systems, which are of great interest in gravitational wave astrophysics fsee lSathvaprakash fc Schut zH2009T). 
Many neutron stars reside in binary systems, e.g. X-ray binaries and pulsar binaries (see Liu et al. 20061 . 20071 : Lorimerl 



2008). Some neutron stars also form pairs, known as double neutron star (DNS) systems, e.g. PSR B1913+1 6. DNS systems 
are expected to be rare. However, more than a dozen DNS systems have already been found (see lLorim er 2008) . The discover y 



of the DNS system B2127+11C, which is probably being ejected from the core of the globular cluster M15 i Prince et al.lll99ll ). 
suggests that dense stellar en vironments could produce DNS systems efficiently. The progenitors of DNS systems are massive 
star binaries. A recent study (|Sana et al.ll2012l ) showed an unexpectedly large fraction of massive stars are in binary systems 
with over 70% of all massive stars transferring material to their companions. Finding such a large fraction of massive star 
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binaries with substantial mass exchange clearly indicates that the evolution of massive binary stars is far from trivial and 
that DNS systems may be more abundant than previously thought. 

Given a substantial number of neutron stars in the central regions of galaxies, some of them would form DNS systems. 
When a DNS system sinks deeply toward the massive central black hole of its host galaxy, it would end up in a close orbit 
around the black hole before the final in-spiral. DNS systems are detectable if they contain a pulsar. However, detecting 
puls ars and binary pu lsars in the centres of galaxies, e.g. the Galactic Centre, is a great technical challenge in observation 
(see lBates etai]|201lft . 

The dynamic s of strongly bound gravitational systems often harbor interesting complex phenomena. Studies (e.g. ISinghl 
120051 ; ISingbl l2008h have shown that a fast spinning neutron star orbiting around a black hole will exhibit complex spin 
precession due to various relativistic couplings. For pulsars, t he spin precession of n e utron stars orbi ting around black holes 
will be manifested in the variation of their pulse emission (|Wex fc Kopeikinlll999l ; iLiu et all [20121 ). Pulsars are therefore 
useful experimental probes of relativistic spin-orbit interactions. Hierarchical three-body systems, such as compact binarie s 
orbiting around massive black holes, are known to show complex dynamical behaviours (see e.g. lAntonini fc Peretsll2012h . 
Tightly bound DNS systems are a special subclass of such hierarchical three-body systems. Like pulsars, DNS systems are 
also spinning objects, and hence they would experience similar spin-orbit interactions as they revolve around a black hole. 
However, DNS systems are not point gyroscopes, unlike pulsars, and the subjection of a DNS system to tidal interactions in 
the gravitational field of the black hole gives rise to richer dynamical behaviours. 

In this study, we explore the complex orbital dynamics in the hierarchical three-body system consisting of a DNS system 
orbiting around a massive black hole. We consider the tightly bound DNS systems, which have orbital periods much shorter 
than the periods with which they revolve around the black holes. The focus of this work is on the complex dynamics arising 
from the interaction between the orbit of the two neutron stars and the orbit of the DNS system around the massive black 
hole. We will leave the more complicated interaction between the orbit of the two neutron stars and the spin of the black 
hole to a future study. We organize the paper as follows. We first formulate the dynamical interaction between the internal 
and external orbits of the DNS system revolving around a black hole and derive the corresponding effective potentials (§2). 
We recast the formulation in terms of Euler angles for the orientation of the (internal) orbital plane of the neutron stars in 
the DNS system and construct the Lagrangian, complete with terms corresponding to Coriolis and centrifugal forces for the 
orbiting reference frame. Then we derive the equations of motion that govern the orbital dynamics of the system (§3). We use 
two special cases to illustrate the precession and nutation behaviours, and finally we present a full numerical solution (§4). 
The physical and astrophysical implications are discussed. 



2 INTERACTIONS AND EFFECTIVE POTENTIALS 

The configuration of the hierarchical DNS and black hole system is shown in Fig. [0 We treat the black hole and the neutron 
stars as point masses. As neutron stars have very similar masses, around 1.5 M© " (|Lattimerll201lh . we consider that the two 
neutron stars have equal masses m. The mean orbital separation between the neutron stars is a, and the distance of the 
neutron star to their centre of mass is therefore a/2. The DNS system is orbiting about a massive black hole with a mass M 
substantially larger than the that of the neutron star mass, i.e. M m. The orbital radius of the DNS system is R. Moreover, 
R 3> a. As the primary aim of this study is to demonstrate that complex behaviours can arise from orbital coupling in the 
DNS systems even in very simple configurations, we avoid unnecessary complications. We consider Schwarzschild spacetime 
for the black hole and zero orbital eccentricity for the DNS system. 

The orbital angular frequency of the neutron stars in the DNS system is given by Q nB = J Gm/2a 3 , and the orbital angular 
frequency of DNS system around the black hole is given by Q. bh = s/GM/R 3 . The condition of R/a 3> (2M/m) 1//3 implies 
that f2 ns ^> fibh- This condition is easily satisfied with the systems of our interest, so the DNS system can be considered 
as tightly bound. The fast orbital revolution of the neutron stars around each other effectively generates a looped mass 
current. On the time scales t 3> 27r/fibh, the continuous approximation is applicable and the looped mass current effectively 
makes the tightly bound DNS system into a relatively rigid rotating mass ring. It is important to note that a two-body 
DNS system is not identical a rigid rotating mass ring, as the orbit of the former is maintained by gravity and the shape of 
the ring is determined by the constraint force. However, in the Lagrangian formulation for the dynamical analysis of a rigid 
ring in an external field, consideration of constraint forces within the ring need not be included explicitly. In essence, the 
distance a/2 of each diff erential mass element from the centre of the ring forms a holonomic constraint (see [Goldstein ll95Cll ; 
iFowles fc Cassidavll2005l), with the re sult that any displacement from this rigid configuration is an ignorable coordinate in the 



Lagrangian (|Fowles fc Cassidavll2"005h . Thus, the justification of the the rotating mass ring representation for the DNS systems 



gives us tremendous simplification in the mathematics when carrying out the analysis of the dynamical coupling between the 
internal orbit of the neutron stars in the DNS system and its external orbit of the DNS system around the massive black hole. 



2.1 Tidal interaction 

We start the derivation of the effective potential by considering the tidal force on the DNS system. The tidal force between 
freely falling objects moving along two adjacen t geodesies with separation r) m ay be derived from the Riemann curvature 
tensor in the local, orthonormal, inertial frame (|Misner et al.ll 19731 ; [Rindlerll200d ) : 
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Figure 1. A schematic illustration of the configuration of a tightly bound double neutron star (DNS) system orbiting around a massive 
black hole (not to scale). The orbital separation of the two neutron stars is a and the radius of the orbit of the DNS system around 
the black hole is R. The orbit of the neutron stars in the DNS system and the orbit of the DNS system around the black hole are not 
necessarily co-planar. 
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where dr denotes the unit of proper time at that spacetime location and D is the covariant derivative operator. For the 
Schwarzschild metric, 
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(see also lFang fe Ru ffini 1983|). The tidal force stretches objec ts in the radial direction and compresses them within the plane 
normal to the radial. As derived in lOhanian fe Ruffinil <|l994h for suitably small as in the case of the DNS and black hole 
system considered here, where a <C R, Eq. [T]may be expressed in the Schwarzschild time coordinate t. Hence, Eqs. [2] through [4] 
may be recast in terms of an effective potential, seen by a distant observer, acting relative to the DNS system's centre of mass: 



$ = _ 



1 GM 
1 GM 



2{rf 



(5) 



since (jf) 2 + (r) e y+ (rf) = (a/2) 2 . 

The equipotential surfaces of $ are hyperboloids, with a saddle point at the centre of mass of the DNS system (the mass 
ring), as illustrated in Fig. [21 This potential is deriv ed in the weak-field, slow-mo tion approximation. It is a special case of the 
gravitational quadrupole potential as that shown in lBarker fe O'Connelll (|l975al ). Here we have shown that for a DNS system 
orbiting a black hole the potential in the mass ring representation may be simply derived as following immediately from the 
equations for the tidal force. 

The potential energy Vtidai of the mass ring as a result of its orientation within the tidal field is given by integrating $ 
over each mass element in the mass ring, i.e. 



(6) 



where A = 2m /na. Let the plane of the ring be defined by normal unit vector n, with zenith and azimuth coordinates 9 and 
<j>. Parameterizing the mass elements of the ring in terms of an angular element d\ — 2 dl /a permits the expression of rj' for 
a particular element via coordinate rotation of the vector (cosx, sm X) 0) T , which describes a ring in the x-y plane, into the 
plane of the ring: 
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Figure 2. Potential contours, illustrating slices through the effective potential <!> of the Schwarzschild tidal force. Potential increases from 
dark to light shading. Left: Hyperbolic contours, of a saddle-like surface, in a plane containing the radial vector. Centre: Circular contours, 
of a paraboloidal surface, in the plane normal to the radial vector. Right: Hyperboloidal equipotential surfaces in three dimensions. 
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Here the azimuth angle is expressed in co-rotating coordinates, cj>' = cj) — ilt, where fi (= fibh = uGM / R 3 ) is the angular 
speed of the orbit in coordinate time. Hence, (/>' is the angle between the radial vector from the central gravitating body (the 
massive black hole) to the ring (the DNS system), R, and the projection of n into the plane defined by R and RxL, where 
L is the orbital angular momentum. That is, when L is oriented along the z direction, so that the orbit about the massive 
body is in the x-y plane, <f>' is the angle between R and the projection of n into the x-y plane. Integrating yields the total 
orientational potential energy of the ring: 

1/ 1 GMm ( a \ 2 W Tl ->( 2„ 2 ,/ 2 . . 2 ,/ . 2 a . ,/ . \" 

vtidal = ^3 — y — J / d-x 1 — 3(^cos 8 cos <p cos x + sm 4> sm X — ^ cos 8 cos (p smq> cosxsmxj 

GMm ( a \2 f 3 r 2 . 2 



fit) + sin 2 (0 - fit)] 



(10) 



Note that Eq. 1101 which we have derived completely classicall y, may be compared with the orientational potential energy 
from the Hamiltonian for the g r avitati onal t wo-body problem in [Barker fc O'ConneiH (|l975bl ') discussed in a semi-classical 
context in iBarker fc O'Connelll (|l975al ) and IChan fc O'Connelll (| 19771 ), which gives the potential energy for a mass with 
quadrupole moment J 2 . In terms of our notation, their expression is 



Kidal = 



GJ 2 M (2m) 
2R? 



3 (n ■ R) 



- 1 



(11) 



For the thin ring, Ji = (1/2) {a/2) 2 , so with R = R (xcos fif + y sin Qt) and n = x sin 8 cos <j> + y sin 8 sin (j> + z cos 0, n- R/_R = 
sin 8 cos <f> cos f2i + sin 8 sin sin Q,t = sin 8 [cos cos (—Qt) — sin sin (—fit)] = sin 8 cos (4> — Oi). Hence, Eg. 1111 becomes 

1 GMm 

2 R 

GMm / a \ 2 f 



f^-) 2 [3sin 2 6'cos 2 (( ? !)-nf) - ll 
V 2R J 



fit) + sin (cf> - ttt)] 



R \2R) 1 2 

identically the result obtained through the effective potential formalism in Eq. 1101 



(12) 



2.2 Orbital coupling 



We next derive the potential due to orbital interaction. The transport of vectors, such as angular momentum, along paths 
in curved spacetime gives rise to precessional effects. Within the Schwarzschild metric, the transport of angular momentum 
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vectors along circular orbits is subject to deSitter precession jde Sittedll916h . That is, the orientation of the vector will precess 
with angular veloci ty proportional to rxv, where r and v are the position and velocity of the test body carrying the angular 
momentum vector ijMisner et al.lll973l : iRindlerlliool l . This implies that when the DNS system revolves around the black hole, 
the orbital angular momentum vector of the neutron stars about each other will gradually change its orientation. 

The deSitter precession can be derived geodesically, i.e. directly finding the deflection resulting from parallel transport, or 
equivalently, by converting the Schwarzschild metric to a rotating coordinate system. However, if a Lagrangian analysis within 
the Schwarzschild metric is required, it is more calculationally useful and conducive to physical insight to express the coupling 
between the angular momentum of the spinning body about its centre of mass (S), and the orbital angular momentum (L) as 
a potential energy, Vsl- Obviously, here S corresponds to the summed orbital angular momentum of the neutron stars in the 
centre of mass frame of the DNS, and L corresponds to the angular momentum of the DNS system's orbit around the black 
hole. 

We now show that the general relativistic orbital coupling between the DNS system and the black hole can be derived 
using the gravito-electromagnetic (GEM) formalism. Note that the results are also applicable for spin-orbit coupling in binary 
systems of arbitrar y mass ratios. The GEM equatio ns constitute a next-to-leading order approximation to the Einstein field 
equations (see e.g. Misner et al.lll973l ; |Punslvll200ll ). They closely resemble the form of Maxwell's equations. In essence, the 



derivation of the gravitational spi n-orbit interac tion energy is analogous to the derivation of the spin-orbit energy of the 
electron in the hydrogen atom (see I Jackson! 1 1 99^ ) . 
The GEM equations read 

V ■ g = 4^Gp g ; (13) 

V -| = 0; (14) 

Vxg = 44 (V) i (15) 



c dt V 2 

_ h 4nG . Ids . 

Vx 2 = — J ^ + ci' < 18 > 

where g is the negative of the usual gravitational field strength, h is the gravito- magne tic field, p g is the mass density, and 
j g is the mass current density (|Punslvll2o"oil ; iRuggiero fc Tartagliall2002l ; |Mashhoon1l2007l '). The GEM Lorentz force law for an 
object of mass m is 

F = -m (g + 2-xh) . (17) 

In the above equations the additional factors of 2, compared to their electromagnetic counterparts, are caused by the second- 
rank tensor structure of the Einstein field equations, in contrast to the first-rank tensor formalism of classical electromagnetism. 
The sign differences are the result of the sign difference between the gravitational and electromagnetic forces among analogous 
configurations of particles (e.g. two positive charges repel one another, but two positive gravitational "charges" attract each 
other). 

For clarity, let 1 and 2 denote the non-spinning and spinning bodies (the massive black hole and the orbiting DNS 
system), respectively. Thus, mi = M and 7712 = 2m (see §2). Let and v; denote the positions and velocities of bodies 1 and 
2, measured in the centre of momentum system. In the reference frame of body 2, the spin-orbit potential energy is given by 

^ ) =-2m 2 .h( 2 ) = -2m 2 .(h«+2^xg( 1 ») , (18) 

where m 2 is the gravitational analogue of the magnetic dipole moment for body 2, v = v 2 — vi is the speed of body 2 relative 
to body 1, and h' 2 ' is the gravito-magnetic field in the frame of body 2, which transforms into fields g' 1 ' and h' 1 ' in the 
frame of body 1. In the rest frame of body 2, a mass current is observed as a result of the apparent motion of body 1 about 
body 2, analogous to the apparent orbit of the proton in the electron's rest frame in the electromagnetic spin-orbit problem. 
Just as the proton's apparent current gives rise to a B field in the electron's frame, the apparent motion of body 1 gives rise 
to nonzero YS 2 ^ . The additional factor of 2 at the beginning of Eq. [TH] is due to the doubling of the h field contribution in 
Eq. 1171 the frame transformation also follows from Eq. 1171 Further, h' 1 ' = 0, since the gravito-magnetic component of the 
field vanishes after h' 2 ^ is transformed to the reference frame of body 1. That is, the gravito-magnetic field created by body 
2 in the frame of body 1 is not relevant, since we are interested only in transforming the fields seen in the frame of body 2. 
Equivalently, 

Vg=-&. e ™, (19) 
where 



pf =4m 2 x- (20) 



v 



c 



is the effective gravito-electric dipole moment created by the motion of m 2 in the frame of body 1. We not e that Eq. [20l is 
directly analogous to the electric dipole p = v x m created by a magnetic dipole m mov ing at velocity v |Fisherlll97lf ); a 
consequence of special relativity, this effect was first predicted bv lEinstein fc Laubl (| 19081 '). In this case, g 1 - 1} = Gmiir 3 R, 
where R = r 2 — ri is the position of body 2 relative to body 1. Incorporating Eq. [20]into Eq. 1191 yields 
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( 2) AGm 1 R 4Gmi 
V sl = —Jp ■ ( ma x v ) = TTrT m 2 ■ (R x v) . (21) 



A familiar result from electromagnetism is the gyromagnetic ratio Y connecting the magnetic moment and angular 
momentum of a body such that m = TS; for a body of mass m and charge q, V = q/2m if q/m is uniformly distributed. The 
gravitational analogue of V can be viewed as being generated by mass currents within the body, due to S. Dividing body 2 
into rings coaxial with S2, its angular momentum about its centre of mass, results in F = 1/2, that is, 

m 2 = is 2 , (22) 

in direct analogy to the classical electromagnetic result (i.e. an object's gravitational charge is identically its mass). Incorpo- 
rating this result into Eq. 1211 yields 



We perform a frame transformation, changing from the reference frame of body 2 to the centre of momentum frame of 
the complete two-body system. In order to evaluate Vsl in the centre of momentum frame, the Thomas precession must be 
taken into account. Expressing this effect as a potential energy, we obtain 

V$ = S 2 • ^ T ,2 , (24) 

where o^x,2 = (&2 x V2) /2c 2 = —Gmi (R x V2) /2c 2 R 3 , with a.2 = -Gmi_R" 3 R being the acceleration of body 2 in the 
centre of momentum frame. This Thomas precession is a special relativistic consequence of frame transformation. It is a direct 
analogue of that obtained in cha nging from the frame of the electron to the centre of momentum frame of the hydrogen atom 
in the electrodynamical problem (|jacksonll 19991 ). In the GEM formulation effects of spacetime curvature are incorporated into 
the effective fields of the GEM equations. By substitution, 

^i?=-^S 2 .(Rxv 2 ) . (25) 

The orbital angular momentum of the system is L = R x P and, to first order, P = /_tv, where /1 = mim.2/ (mi + m.2) is 
the reduced mass. Hence, Eg. 1231 may be re-expressed as 



flC^H^ C z Hr \ 7712 



By definition of the centre of momentum system, mivi = — 7712V2. To leading order, L = m2R x V2. Thus, the Thomas 
contribution, Eq. 1251 may be written as 



^=-2^( S2 ' L >- < 27 > 
The potential energy of the gravitational spin-orbit coupling in the centre of momentum frame is therefore 

+ + <■.■!.>■ W) 

For the case of the DNS system of mass 2m, with angular momentum S about its centre of mass and orbiting a black hole 
with mass M S> 2m, the spin-orbit energy takes the form: 
o (-1 \,r 

V " = 2(2^3 (S-L). (29) 

Note that this result has a one-to-one correspondence with that of IChan fc O'Connelll (|l977l ). despite the fact that it was 
derived from a completely different physical setup, with a finite spinning mass (the DNS system) as opposed to a spinning 
test point particle. 



3 ORBITAL DYNAMICS 

3.1 Lagrangian in the orbiting reference frame 

We choose the location of the black hole as the origin of the observer's coordinates. Let S denote the distant inertial observer's 
frame, and S' denote the frame rotating with angular velocity f2 = i\jGM/ R 3 , with origin at the centre of the DNS system, 
measured in S as R. Let r and r' represent positions measured in S and S' , respectively. The Lagrangian of a particle of mass 
m, expressed in the coordinates of frame S' and coordinate time t, is 

C = imlr' + fix (r' + R)| 2 -m$(r') + ^^-F Si 

= \m{\v'\ 2 + 2fi- [(r'+R) x r'] + [fiV - (fi ■ r) 2 ] } - m* (r') + - V Si , (30) 



Orbital dynamics of a double neutron star system around a black hole 7 



where m$ (r') is the tidal potential energy of the particle and Vsl is the GEM spin-orbit coupling. The dot operator denotes 
differen tiation with respect to t, i.e. ' = d/dt. For discussions on the Lagrangian formalism in non-inertial reference frames, 
see e.g. iDallen fc Neuenschwanderl (|201lh . 

Here we do not consider explicitly the gravitational interaction between the two neutron stars. As the DNS system is 
represented by a spinning mass ring, the gravitational interaction between two neutron stars is analogous to the constraint 
forces within the ring. As discussed in §1 the distance a/2 of each mass element from the centre of the ring forms a holonomic 
constraint. A displacement from the configuration is an ignorable coordinate in the Lagrangian and therefore there is no need 
to include it in the Lagrangian for the dynamical analysis of a ring. Nevertheless, the inclusion of such displacements in the 
Lagrangian would allow one to calculate the tidal strain within the ring and the degree to which the ring flexes for internal 
interaction, which is an entirely different problem from that consider ed here. F o r mor e detailed discussions on t he role of 
holonomic constraints in the Lagrangian formulation of mechanics, see iGoldsteinl ||l950h : lFowles fc Cassidavl (|2005l ). 

For the ring as a whole, of mass 2m, the Lagrangian may be found by integrating Eq. [30] over the mass elements dm: 

^ 2GMm , f dm 



R 
Now, 



+ f ^{\r\ 2 + 2n-[(v'+-R)xv'] + [n 2 r 2 -(n-r) 2 ]}-VsL-f dm$(r') . (31) 



ing ^ 



dm 'rf =T rot , (32) 



the (rotational) kinetic energy of the ring in S'. Further, 

dm$(r') = Kidai , (33) 



ring 



where Vtidal is the orientational potential energy created by tidal forces. A further simplification of Eq. [31] is achieved by 
expanding SI ■ [(r' + R) x r'] = SI ■ [(r' x r') + (R x r')]. Since SI ■ (R X r') = f ' • (SI X R) and SI x R does not depend on r', 
the integral 



ring 



(R x r) = (SI x R) ■ * dm r = , 

" ring 



dm SI ■ (R x r) = (SI x R) • * dm r = , (34) 



by symmetry. Thus, the Lagrangian becomes 

C = Trot - Kidai - Vsl + + / ^ {2SI ■ (r' x r') + [ftV - (SI ■ r) 2 ] } . (35) 

iring ^ 

The physics is more transparent if the problem is expressed in terms of Euler angles. Let the 123 coordinate system be 
the body coordinate system moving with the orientation of the ring, with the 3-axis defined as the symmetry axis. Define 
the S" coordinate system as follows. Let the x"-axis correspond to the line of nodes, that is, the intersection of the 1-2 plane 
(plane of the ring) with the x-y (equivalently, x'-y') plane. Let z" be aligned with the 3-axis. Finally, let a' , /3', and 7' denote 
the three Euler angles for the ring in frame S' . That is, let a' be the angle between the z- (equivalently, z'-) and z"-axes, ft' 
the angle between the x'- and s"-axes, and 7' the angle between the x"-axis and the 1-axis, describing rotation of the ring 
about its symmetry axis. In terms of n, the normal vector to the plane of the ring, with zenith and azimuth coordinates 6 and 
<f>, a' = 9 and /?' = cj> — Qt + tt/2. In terms of the Euler angles defined in the inertial frame 5, where ft is the angle between 
the line of nodes and the x-axis, a' = a, ft' = ft — Sit, and 7' = 7. Therefore, hereafter in the analysis, primes will be dropped 
from a and 7. 

The angular velocity of the spinning mass ring, which is also the angular velocity of the neutron stars, in S' may then be 
expressed as u>' — a + 0' + 7. In terms of the 123 coordinate system, 

Lj[ = p sin q sin 7 + d cos 7 ; 
LJ2 = p sin q cos 7 — d sin 7 ; 

Lj' :j = p cos a + 7 . (36) 

The kinetic energy T Iot is given in 5' by (1/2) $^ =1 hu' 2 , where h = I 2 = (1/2) (2m) (a/2) 2 = ma 2 /4 and J 3 = (2m) (a/2) 2 = 
ma 2 /2. Thus, 

2 f /i/2 2 2 22 "/ " /2 2 2 2 2 ' f 

Trot = —ma 1/3 sin a sin 7 + d cos 7 + 2a/3 sin a sin 7 cos 7 + p sin a cos 7 + a sin 7 — 2d/3 sin a sin 7 cos 7 
8 V 

+2/3' 2 cos 2 a + 27 2 + 4/3' 7 cos 

= ima 2 [/3' 2 (1 + cos 2 a) + a 2 + 2 7 2 + 4/3' 7 cos a] . (37) 
The non-inertial terms in the Lagrangian can be divided into a Coriolis term 
Tcor = f dmSl- (r X r) (38) 
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and a centrifugal term 

T ccn = ~ I dm [fiV - (fl ■ r) 2 ] . (39) 

J ring 

The Coriolis term may be simplified by expressing r' as x r'. We may rewrite ft ■ (r' x r') as 

fl • [r' x (w' x r')] = fl ■ [J (r' ■ r') - r' (r' ■ w')] = fl ■ J r 12 - (fl ■ r') («' ■ r') . (40) 

For the rigid mass ring, r' = a/2 is a constant. In terms of the S" coordinate system, fl = fl (y" sin a + z" cos a) 
and a/ = x"d + y"/3'sina: + 2" f/3' cos a + 7^. For a mass element dm along the ring, specified by the Euler an- 
gle 7, r' = r' (x" cos 7 + y" sin 7) . Thus, fl • r' = fl r' sin a sin 7, a/ • r' = r' (a cos 7 + p" sin a sin 7^ , and fl • w' = 
f2 (ji' sin 2 a + $' cos 2 a + 7 cos = SI f /3' + 7 cos . Therefore, expressing dm as 2m 07/271-, 

Tbor = 2m f — ") fl ■ w' — - — ^ — ^ f2 ^ c?7 (sin a sin 7) ^d cos 7 + sin a sin 7^ 



2/ 2tt 



1 2 o 
— ma iZ 



j8' ( 1 — — sin 2 a ) + 7 cos a 



(41) 



Further, the centrifugal term may be simplified by expansion. Since r = R + r', we have 

SlV - (fl • r) 2 = fl 2 (R + r) ■ (R + r') - [fl • (R + r')] 2 

= fl 2 (R 2 + r' 2 + 2R ■ r') - (fl • R) 2 - (fl ■ r') 2 - 2 (fl ■ R) (fl • r') . (42) 

Now, fl _L R and / ring dm R ■ r' = R ■ / ri dm r' = 0, by symmetry. Thus, with fl ■ r' = Sir' sin a sin 7, as shown previously, 
the centrifugal contribution to the Lagrangian is 

= ^ 4> dm |flV* - {fl-r')"\ + mfl 2 R 2 



1 I dm \fl 2 r' 2 - (fl ■ r') 2 ] 

2 /ring ^ J 

i (*», (I)' tf + rfrf - i (£) (!)' tf *"„jf <W 7 



2 

= -ma 2 !! 2 ^1 — — sin 2 c^j + mfl 2 R 2 . (43) 

Upon inspection of Eq. 1431 one can see that in fact the centrifugal term can be interpreted as a potential energy in the 
Lagrangian, with V cen = -T cen . 

The tidal orientational potential energy was derived in Eq.[l0]in terms of the direction of n. It becomes, in terms of the 
Euler angles, 

V tid ai = — r— (al) |l- - [cos 2 asm 2 /3' + cos 2 /?'] J . (44) 

Finally, the relativistic spin-orbit potential energy must be expressed in terms of the Euler angles. The orbital angular 
momentum is L = 2mR 2 fl and the angular momentum of the ring is S = Iw, where I is the ring's moment of inertia tensor 
and oj is the angular velocity of the ring defined in S. With denoting the orthonormal basis vectors of the 123 coordinate 
system, oj = ei (j) sin a sin 7 + d cos 7^ + §2 L8 sin a cos 7 — d sin 7^ + §3 N3 cos a + 7^ . Note that here we use /j, not /?', since 
it is the spin observed in the inertial system that determines Vsl- Working in the 123 coordinate system, we have: 

S = ~ (2m) ^ ^ ^/3sin«sin7 + d 0087^ ei + ^/3sinacos7 — asin7j e 2 + 2 cos a + 7^ e 3 j (45) 

and 

fl = fl (ei sin a sin 7 + 62 sin a cos 7 + 63 cos a) , (46) 
whence 



SL = l -(2m) 2 (^) 2 R 2 fl 



* 2 2 * 2 2 "2 

p sin a sin 7 + d sin a sin 7 cos 7 + p sin a cos 7 — d sin a sin 7 cos 7 + 2/3 cos a + 27 cos a 



^m 2 a 2 R 2 fl [/3(l + cos 2 a) +27 cos a . 



Thus, 



3 GM 
F5L = 2(2m) C 2 J R3 ( S ' L ) 



3 GMm a z fl 
8 _ R IF 



\ (f}' + fl\ ( 1 + cos 2 a) + 2 7 cos a] . (48) 
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Finally, the complete Lagrangian for a DNS system of mass 2m and diameter a in circular orbit at distance fl> a from 
a black hole of mass M S> 2m may be written: 

C = \ma 2 [/3' 2 (1 + cos 2 a) + a 2 + 2 7 2 + 4/3' 7 cos a] [T TO t] 

+ ±ma 2 Q, \p' (1 - | sin 2 a) + 7 cos a] [T C or] 

+ima 2 tf (1 - i sin 2 a) [Ken] 

-^{mf [l-|(cOS 2 QSin 2 /3'+COS 2 /?')] [Kidal] 

_|Gm™^ [(/?' + + cos 2 q) +2 7 cosa] [Vsl] 

+mQ. 2 R 2 + , iGMrn L _ [constants] (49) 



In terms of f2 = yjGM/R 3 where possible, we may express the Lagrangian as 
C = \ma 2 [/3' 2 (1 + cos 2 a) + a 2 + 2 7 2 + 4/3' 7 cos a] [T rot ] 
+ ima 2 fi \/3' (1-1 sin 2 a) + 7 cos a] [T Cor ] 

+ ima 2 n 2 (1- isin 2 a) [Ken] 
-|ma 2 fi 2 [1 - § (cos 2 asin 2 ft + cos 2 /?')] [Kidai] 

-fmfi 3 ^ [(/3' (l + cos 2 a) +2 7 cosq] [Vsl] 

+3mQ 2 R 2 . [constants] (50) 

Grouping like terms, applying trigonometric identities, and rewriting Eq. [50] in terms of /3 = + ^i, that is, in the Euler 
angles of the inertial frame, we obtain 

£ = ima 2 [/? 2 (1 + cos 2 a) + a 2 + 2 7 2 + 4$j cos a] [T rot , T C or, Ken] 

-±ma 2 fi 2 {l- § [cos 2 a sin 2 (/3 - fit) + cos 2 (/3 - fit)] } [K tidal ] 
_| m f2 3 ^ [/3(1+cos 2 q) +2 7 cosa] [K Si ] 

+3mQ 2 R 2 . [constants] (51) 

Trot, Tcor, and Ken combine to yield the rotational kinetic energy seen by an observer in the inertial frame. Through this 
analysis, we have shown explicitly that the Lagrangian can be correctly formulated in the inertial frame, provided that the 
tidal potential seen in the non-inertial frame undergoes time-dependent coordinate rotation at rate f2. Coriolis and centrifugal 
forces do not otherwise affect the results. 



3.2 Equations of motion 

The Euler-Lagrange equations governing the dynamics of the DNS system are 

±(K\-K ( 52 ) 

where ip = a, /3, 7 and, as previously, the dot operator signifies differentiation with respect to coordinate time, ' = d/dt, 
that is, the time measured by a distant inertial observer. Proper time, experienced in the reference frame of a particle, is 
defined differentially as cdr = yj dx^dx^ = yj g^dx^dx" , where <? M „ is the spacetime metric tensor. The Schwarzschild metric, 
expressed in spherical coordinates, is 

_ Ti\ c * dt 2 - (l - -) ~* dr 2 - r 2 d9 2 - r 2 sin 2 8d0 2 , (53) 

where r B = 2GM/c 2 is the Schwarzschild radius. For particles at rest at some radius R in the Schwarzschild metric (dr = 
dO = dcj> — 0), dr = dt\J\ — 2GM/c 2 R. For particles in circular orbit at that radius, dr = d6 = and 6 — tv/2, but d<j> = Qdt. 
Here, £7 is the u sual Keplerian ang ular velocity, y/GM/R 3 , which is the exact general relativistic result in coordinate time, 
regardless of R jMisner et al.|[l973r ). Thus, for the orbiting particle, 

1 2GM \ 2 , 2 GM 2 A 3GM\ 2 , +2 ,_,v 
It follows that 
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Taking the relevant derivatives of Eq. [51] and rearranging, we obtain the equations of motion for the DNS system: 



= a+ ^/3 2 -3fi 3 — /3J sin a cos a + ^2/3-3fi 3 -^-J 7 sin a + 3fi 2 sin 2 (/3 - fit) sin a cos a ; (56) 
( R 2 \ 

= '$ (l + cos 2 a) — f 2/3 — 3fi 3 — — J a sin a cos a — 2ccy sin a + 2j cos a + 3fi 2 sin a sin (/3 — fit) cos (/3 — fit) ; (57) 

= 27 + 2/3 cos a- ^2/3 -3fi 3 ^-^ a sin a . (58) 
Defining a new parameter £ = 2/3 — 3fi 3 -R 2 c -2 , Eqs. 1561 through 1581 can be expressed more compactly: 

= a + ( £ — f3^ f3 sin a cos a + £7 sin a + 3fi 2 sin 2 (/3 — fit) sin a cos a ; (59) 

= /3 (l + cos 2 a) — £a sin a cos a — 2A7 sin a + 27 cos a + 3fi 2 sin 2 a sin (/3 — fit) cos (/3 — fit) ; (60) 

= 27 + 2/3 cos a — £d sin a . (61) 



4 RESULTS AND DISCUSSIONS 

4.1 Demonstration of restricted oscillations 



To understand some of the dynamical properties of this system, it is instructive to consider the special cases in which the 
tidal effects are singled out. If we let the mass ring (the DNS system) be held at fixed location in space by a force, the time 
dependence of Eq. [51] is suppressed. Without loss of generality, we fix the centre of mass of the ring to lie along the positive 
i-axis, so that the orientational potential energy, now time-independent, may be expressed in terms of the coordinates of n: 



Kid 



GMm 



(—)' 
R \2RJ 



> cos 2 d> + sin 2 1 



(62) 



The behaviour of these special, restricted cases will be useful in understanding the more complex nature of solutions to the 
general case of the orbiting, unconstrained ring. The Lagrangian for this restricted system is 



stat — 



1 2 

—ma 

8 



+ (1 + cos 



GMm / a \ 
R \2r) 



I- 3 - 
2 



' cos ffl + sin 



(63) 



which implies that the oscillations have no explicit dependence on the orbital parameters of the neutron stars in the DNS 
system, in particular, the orbital separation a. 

We consider two cases to demonstrate the restricted oscillations arising from such situations. In the first case, we let <j) 
be fixed and let 9 be free; in the second case let be fixed and let <f> be free. In the calculations, the system parameters are 
M — 1.5 x 10 7 Mq (black hole), m = 1.5 Mq (neutron star), R = 25 r B , where the Schwarzschild radius of the black hole 
r s = 2GM/c 2 = 4.4 x 10 12 cm, and fi = y/ GM / R 3 = 3.8 x 1(T 5 rads" 1 . 



4-1.1 Case 1: fixed <j>, free 9 
Lagrange's equation yields 

9 + 3fi 2 cos 2 (j> sin 9 cos 9 = . (64) 

By inspection of Eq. 1641 one can see that 9 has stable equilibria at integer multiples of n and unstable equilibria at odd 
half-integer multiples of ir. It is easy to visualize in the <f> = case: tidal forces will try to pull the mass ring into the plane 
defined by R and R x L, that is, the x-y plane. For small 9, sin 9 cos 9 sa 9 and oscillation about stable equilibrium will occur 
with angular frequency 



u>0 ,0 



fi\/3|cos0| , (65) 



that is, 9 = 9q cosoje,ot- 

This estimate may be improved by the method of successive approximations, using the next term in the Taylor expansion, 
cos 9 sin 9 m 9 — (2# 2 /3) , and the result obtained in the first approximation. This process results in a slight lengthening of the 
period of oscillation and the appearance of higher harmonics. At next-to-leading order, 



ujg = u g ,o\j 1 - ^#0 (66) 
and 



= cos ui e t + (276»o - 48) 6»o cos 3w e t 



(67) 
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Figure 3. Phase space and time-dependent behaviour in the zenith angle 6 for oscillation of the fixed, 0-restricted, non-rotating ring. 
For small oscillations, the behaviour is simply harmonic, while for oscillations with amplitude approaching it/2, lingering behaviour 
occurs as a result of the unstable equilibrium. In the calculation we use M = 1.5 X 10 7 Mq (black hole), m = 1.5 Mg (neutron stars), 
R = 25r s , where r s = 2GM/c 2 = 4.4 X 10 12 cm, and ft = ^/gM/W = 3.8 X 1CT 5 rads" 1 . The units of 6 arc rads" 1 . 



In phase space (6 — 6), when 60 is small, the path traces an ellipse, as required of simple harmonic oscillation. However, 
as #0 approaches it/2, n tends to hang longer near this point of unstable equilibrium, making the time-dependent behaviour 
of 9 approach a square-like wave. Likewise, the phase plot becomes somewhat pinched and the phase trajectory approaches a 
homoclinic orbit as the initial condition approaches unstable equilibrium (see Fig. [3j) . 



4-1.2 Case II: fixed 6, free 4> 

The application of Lagrange's equation to Eg. 163 1 yields 

(1 + cos 2 6) 4> - 3Q, 2 (1 - cos 2 6) sin 4> cos <j> = . (68) 

As one can see, Eq. [68] possesses stable equilibria at odd half-integer multiples of n and unstable equilibria at integer multiples 
of 7t. As in the case of 6, small oscillations in <j> about the stable equilibria may be treated as simply harmonic, with natural 
frequency 



w *.°=Wtt£4- (69) 

This case, too, may be corrected through successive approximations, with identical results, for <j> near -n/2: 



u$ = «Wl- i<#) ( 70 ) 
and 

4> = — + 4>o cos uj^t + (2700 — 48) 1 cos 3o^t , (71) 

where <j)o = <f> (t — 0) — n/2. The behaviour of cj> in time and phase space is analogous to that shown in Fig. [3] but shifted by 
tt/2. 

Thus, the natural oscillatory frequency of the mass ring due to tidal forces scales with the orbital frequency fl. This 
relationship causes the orbit to function as a driving frequency for the case of the mass ring in orbit about the central mass. 
In the reference frame of the ring, the effective tidal potential of Fig. [2] is spinning at frequency il. However, this resonant 
behaviour is highly dependent upon initial conditions and the spin of the mass ring. 
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Figure 4. Top left: Numerical solution for the normal vector n to the plane of the spinning ring, for 25 orbits about the black hole. Note 
the nutation, upon which smaller processional effects arc superimposed. The parameters used for this simulation are M = 1.5 X 10 7 Mq 
(black hole), m = 1.5 Mq (neutron star), Ft. = 25r s , where r s = 2GM/c 2 = 4.4 X 10 12 cm, and Q. = ^/GMJW = 3.8 X 10~ 5 rads _1 , 
with the initial conditions a (t = 0) = 0.02 rad, /3 (t = 0) = 7r/2, 7 (t = 0) = 0, a (t = 0) = /3 (t = 0) = 0, and 7 (t = 0) = v / 2Gm/o 3 = 
3.8 X 10~ 4 rads _1 , where a = 2 Rq. Bottom left: Magnification of solution for first full orbit about the black hole, illustrating the 
precessional effects within the nutation. Top right and bottom right: The corresponding results from pure Newtonian consideration are 
shown for comparison. They are different from the results obtained from the relativistic derivations. The Square Kilometer Array (SKA) 
is expected to achieve timi ng precision at the level of ~ 100 ns at 1.4 GHz for 10 min integration for millisecond pulsars with normal 
brightness (jLiu et alj|201lh . The differences in the nutation and precession shown for the two cases can be distinguished, if one of the 
neutron stars in the DNS system is a millisecond pulsar. 



4.2 General situation 

Finally, we demonstrate the nutational and precessional behaviour of the general case. A numerical solution of Eqs. [S5] 
through[58]is calculated for the same parameters as previously: M — 1.5 X 10 7 Mq (massive black hole), m = 1.5 Mq (neutron 
star), R = 25 r B , where r B = 2GM/c 2 = 4.4 x 10 12 cm, and Q = ^/GM/R 3 = 3.8 x 1CT 5 rads" 1 , with the initial conditions 
a(t = 0) = 0.02 rad, /3 (t = 0) = tt/2, 7 (t = 0) = 0, a (t = 0) = $ (t = 0) = 0, and 7 (t = 0) = ^2Gm/a? = 3.8 X 10~ 4 rad s _1 , 
where a — 2 Rq . These conditions were chosen to illustrate the effect of a small angular perturbation on the orbiting system; 
7 (t = 0) was chosen to represent the Keplerian angular speed of the analogous DNS system for the mass and size of the ring. 

The numerical solution for the spatial position of n for 25 orbits about the black hole, as seen from the positive z-axis, 
is shown in Fig. 0] Nutational behaviour is clearly evident, with other precessions superimposed. For this case, the orbital 
motion of the ring about the black hole is approximately 3 percent relativistic; that is, dr ~ 0.97dt. However, the importance 
of relativistic orbit-orbit coupling in dictating the dynamics of the system is much higher; assuming 7 > ^, we find for this 
system that 

.2 



VSL 



Kid 



■6^-«1.2. (72) 



Note that precessions and nutations of the internal orbits of the neutron stars in the DNS system revolving around a black hole 
occur even in a pure Newtonian consideration. However, the dynamics in the relativistic and the pure cases are qualitatively 
different, in terms of the precession frequencies and in terms of the interplay between nutation and precession, which can be 
seen in Fig. [4] 

Greater intuition for the effect of the tidal force in governing the behaviour of the DNS system (the spinning mass ring) 
may be obtained by plotting n within the effective potential energy Vtidal + Veen in co-rotating coordinates 9, 4>' . This is done 
in Fig. [5] for three orbits about the black hole. Note that the ring's oscillations are in phase with the frequency of the rotating 
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Figure 5. Phase trajectory of n (solid line), the normal vector to the spinning ring, plotted in co-rotating coordinates 6 and cf>' = ifr — Qt 
for three orbits about the black hole. Contours (grayscale, dashed lines) indicate the effective potential energy Kidai + Veen, with lighter 
tones indicating higher potential energy. 



effective potential, a resonance suggested by our calculations in Eqs. l65l and l69l The ring is driven to low values of 9 at integer 
multiples of n by the potential peaks and swings to higher values of 9 when it is energetically permissible to do so, i.e. at 
half-integer values of n. 

So far, we have not included the gravitational radiation loss in our analyses. For this system, the effects of gravitational 
waves on the dynamics are actually negligible. This can be understo od as follows. The emission of gravitat ional radiation 
would impact the system on timescales of i gra v ~ c?a A /G 3 m 3 or longer (|Misner et al ]| 19731 ; iFang fc Ruffinill983h . However, the 
phenomena of the system operate at frequencies on the order of f2 = \l GM/R 3 or Q.sl ~ GMfl/ Rc 2 . For the astrophysically- 
motivated parameters used above, f^ grav ~ 10 13 and fis^igrav ~ 10 12 . The radiative loss timescales of gravitational radiation 
are therefore much larger than any dynamical timescales that we study. 



4.3 Remarks 



Strongly bound gravitational systems often exhibit intriguing and complex dynamical phenomena. In particular, gravitational 
spin-orbit coupling and the precessional effects it engenders provide an intriguing arena upon which post-Ne wtonian effects 
such as gravito-magnetism come into pla y. Experimentally confirmed to ~ 0.1% precision by lunar-ranging data (jMurphv et al.l 
120071 ). detected in binary-pulsar sy stems (jStairs et al 12004 1 Breton et al ]|2008l ). and recently investigated by the Gravity Probe 



B experiment (|Everitt et al.l 20111). th e spin-orbit interaction suggests deep analogies between gravitation and many other 



physical processes (see e.g. iBlaul I20I1I , for a discussion of experimental tests of the gravitational spin-orbit effect). Across 
many physical systems, the phenomenon of energetic coupling between different components of the angular momentum is 
of paramount importance. More specifically, the coupling between orbital and intrinsic angular momenta is an important 
effect in both quantum and classical systems. In quantum physics, the spin-orbit correction to the Hamiltonian is crucial to 
understanding the fine structure o f atomic spectra , while in optics and condensed matter physics, respectively, examples include 
the c ase of polarization of light (|Niv et a.1.1 l2008l : iBliokh fc Desvatnikovl [20091 ) and the Joseph son current (iDell'Anna et al.l 
120071 ) In classical physics, s pin-orbit coupling occu rs in fields from Maxwellian electrodynamics (|jacksonl|l999l ) to astronomy, 
in the case of tidal locking (|Escribano et al. I2008T ) . The genera l relativistic problem of the spinning point particle is treated 



m the case ot tidal locking (fLscriban o et alJUUUa ) . 1 he genera l relativistic problem ot the spinning point particle is treated 
perturbatively for orbits about a Kerr black hole in lSinghl (|2008l ). A problem related to precession, involving net rotational and 
translationa l displacements produced by cyclic motion of quasi-rigid bodies (analogous to swimming) on curved manifolds, is 
discussed in I Wisdom! (|2003l ). 

For extended nonsp herical bodies orbiting a central mass, tidal forces create additional precessional, nutational, and 
tumbling effects (see e.g. iMashhoon fc SinghlfeoOfj ). In the non-inertial reference frame of the orbiting body, fictitious forces 
must be taken into account to describe the evolution of the orientation. The behaviour of an orbiting body is governed by 
the combination of spin-orbit coupling, tidal forces, and non-inertial effects, and cannot be described by deSitter precession 
alone. This wide and complex range of gravitational phenomena would occur in astrophysical systems such as a DNS system 
infalling into a massive black hole. We have shown that the orbital revolution of the neutron stars in the DNS system is 
dynamically a spinning mass ring of matter, of mass 2m and diameter a, in orbit at radius R about a central black hole with 
mass M. When the neutron stars are revolving around each other, the DNS system constitutes a (nonspherical) gyroscope. 
The fixed radius a/2 in the ring representation of the DNS system forms a holonomic constraint in the Lagrangian analysis. 
As a consequence, internal restoring forces counteracting the tendency of the ring to flex or stretch induced by the tidal 
interaction do not impact the dynamical behaviour of the ring. Thus, the behaviour of a rigid mass ring orbiting around a 
gravitating object models the relevant effects in relativistic astrophysical systems, such as a DNS system that itself is in orbit 
about a supermassive black hole in a hierarchical binary-structured three-body configuration. 

In general, a non-rotating black hole would no t give rise to a n explicit magnetic term under the GEM prescription (for the 
discussion of black hole gravito-magnetism see e.g. |Punslv||200ll ). The hierarchical binary-structured three-body configuration 
consisting of a DNS system orbiting around a non-rotating black hole is a unique astrophysical setting in which GEM effects 
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arise without invoking black hole rotation (Kerr spacetime). The revolution of the DNS system around the central object 
and the revolution of the neutron stars around each other give rise to mass currents, producing magnetic moments and their 
magnetic-type coupling. This situation is analogous to induced magnetism. In the reference frame of the rigid ring, the black 
hole appears as a mass current (similarly to the manner in which the proton appears as a charge current to the electron in 
the H atom in the electromagnetic problem). Hence, in the DNS system's reference frame, a gravito-magnetic field is present. 
When calculating the potential in this reference frame and changing to the centre of mass frame, Thomas precession will 
appear as a consequence. 

Now, the remaining question is: what are the astrophysical consequences? Pulsars are high-precision spinning gyroscopes. 
The y have been considered as useful tools t o probe the spacetime properties of the massive central black hole of our Galaxy 
fsee lWex fc Kopeikinlll999l; IL" iu et al. Il2012l). Moreover, they can also be u sed to investigate various general relativistic spin- 
orbit ^corrpdhigs(|We^ Like pulsars, DNS systems are also high-precision 



spinning gyroscopes. Hence, DNS systems can also be used in experiments on gravity around black holes. In addition, as DNS 
systems are not point masses, they can probe additional phenomena driven by the tidal gravitational field and coupling caused 
by the mass current loop within the DNS systems. Although the chance that a DNS system with a pulsar is currently in 
close orbit with the massive black hole in the Galactic Centre may not be substantial, as future radio telescopes, in particular 
the Square Kilomete r Array (SKA), will be able to search fo r pulsars around central black holes in galaxies beyond the 
Local Group (see e.g. ISmits et al.l l2009: Ri dley fc Lorimerll201Cf ). such hierarchical three-body systems could be found in some 
external galaxies. 

We note that the analyses in this work are not restricted to DNS systems around massive black holes. They are applicable 
also to any hierarchical three-body systems which contain a massive black hole and a tightly bound double compact object sys- 
tem with various combinations of black holes and neutron stars, e.g. stellar-mass black hole pairs (smBH-smBH), stellar- mass 
black hole and neutron star binaries (smBH-NS), and intermediate- mass black hole pairs (imBH-imBH). The central 0.1 pc of 
ga laxies similar to the Galaxy may contain a few thousand stellar mass black holes and several hundred n eutron stars (Table 
1, Hopman fc Alexander! 120061 ). Taking that roughly 10% of neutron stars would be in binary systems |Antonini fc Peretsl 
|2012| ), one may expect about 10 DNS within 0.1 pc from the centres of these galaxies. This is an optimistic estimate, as 
binaries in galactic centres are subject to various 2 -body, 3-body and other non-linear scattering process es, which would lead 
to binary disassociation and evaporation (see e.g. iHopmanl [20091 ; |Peretsll2009l ; lAntonini fc Peretsl |2012| ). Nevertheless, DNS 
and other compact binaries will exhibit complex orbital dynamical properties due to orbit-orbit coupling and some of these 
systems would exist in the central regions of galaxies. Tightly bound double compact object systems are gravitational wave 
sources and their complex orbital precession and nutation will certainly have observational consequences in the properties of 
the gravitational waves from t hese systems. The discovery of pulsars in smBH-NS systems around a massive black hole (see 
iFaucher-Giguere fc Loebll201lh would provide an opportunity for an elegant double-test of general relativity, i.e. comparing 
gravitational wave signatures with pulsar-timing from the system. 



5 CONCLUSIONS 

The behaviour of a DNS system orbiting around a massive black hole is governed in part by coupling between the orbits of the 
neutron stars in the DNS system and the orbit of the DNS system around the black hole. The orbital motion of the neutron 
stars in the DNS systems drives a looped mass current, inducing gravito-magnetism. Because of holonomic constraints, the 
gravitational pull between the neutron stars does not contribute explicitly to the Lagrangian governing the system dynamics, 
so that a tightly bound DNS system is effectively a rigid spinning mass ring. Working in this approximation framework, we 
have constructed the Lagrangian and extracted the equations of motion for the system. 

The dynamical behaviour of this hierarchical three-body system illustrates various general relativistic effects. The tidal po- 
tential energy was derived via a parameterization of the tidal field as an effective potential. The general relativistic spin(orbit)- 
orbit coupling was derived for an arbitrary mass ratio, by application of the GEM equations. The Lagrangian was formulated 
in the orbiting reference frame. It was shown that Coriolis and centrifugal effects collapse into the kinetic energy measured 
in the inertial frame, while the tidal potential undergoes time-dependent coordinate rotation. The Euler angles are found to 
be related though a system of three coupled, nonlinear, time-dependent, transcendental differential equations. Finally, the 
oscillatory behaviour of two restricted cases was examined in detail and a numerical solution of the general case was explored. 

The motion of spinning bodies of finite extent, such as the DNS systems, within a gravitational field contains many 
sources of rich complexity. The system, though seemingly simple, holds many sources of interesting phenomena, among them 
tidal precession and oscillation, spin-dependent nutation, and deSitter precession arising from the spin-orbit interaction. It is 
interesting to note the similarity of the results for both the gravitational quadrupole potential and gravitational spin-orbit 
coupling between the corresponding quantum and classical systems. 
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